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Abstract 

We examine the tripartite entanglement properties of an optical system using interlinked 
interactions, recently studied experimentally in terms of its phase-matching properties by Bondani 
et al. [M. Bondani, A. Allevi, E. Gevinti, A. Agliati, and A. Andreoni, arXiv:quant-ph/0604002 .]. 
We show that the system does produce output modes which are genuinely tripartite entangled 
and that detection of this entanglement depends crucially on the correlation functions which are 
measured, with a three-mode Einstein-Podolsky-Rosen inequality being the most sensitive. 

PACS numbers: 42.50.Dv,03.65.Ud,03.67.Mn,42.65.Lm 
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I. INTRODUCTION 



Entanglement is a property which is central to quantum mechanics and helps to distin- 
guish it from classical mechanics. A vast amount of work has been undertaken on discrete 
variable entanglement, with somewhat less having been performed on the continuous variable 
case. It is the latter which interests us in this work, particularly as regards the entanglement 
of three optical modes. We will focus on an experimentally realised system which links two 
interactions in a combined downconversion and sum frequency generation process Q], 
and examine its utility for the production of states which exhibit full tripartite entangle- 
ment. As far as we are aware, full tripartite entanglement has only been unambiguously 
demonstrated by mixing squeezed vacua with linear optical elements , although other 
methods which create the entanglement using an actual nonlinear interaction are under 
.nvesfgat™, Using both cascaded and _ t ^ processes QflflQ. 

The definition of tripartite entanglement for three-mode systems is a little more subtle 
than that for bipartite entanglement, with different classes of entanglement having been 
defined, depending on how the system density matrix may be partitioned ^. The classifi- 
cations range from fully inseparable, which means that the density matrix is not separable 
for any grouping of the modes, to fully separable, where the three modes are not entan- 
gled in any way. For the fully inseparable case, van Loock and Furusawa Q], who call 
this genuine tripartite entanglement, have derived inequalities which are easily applicable 
to continuous variable processes. More recently, Olsen et al. [ll]] have defined three-mode 
Einstein- Podolsky- Rosen (EPR) ll| type criteria, which also provide sufficient, but not nec- 
essary, conditions for the demonstration of genuine tripartite entanglement. In this article 
we will begin by reviewing the definitions of these entanglement criteria and then apply 
them to the outputs of the Bondani scheme to quantify entanglement correlations which 
may in principle be measured experimentally. 



II. CRITERIA FOR TRIPARTITE ENTANGLEMENT 

We begin by giving the optical quadrature definitions we will use in our analysis, as the 
exact form of the inequalities will depend on these. For three modes described by the bosonic 
annihilation operators a,j, where j = 1,2,3, we define quadrature operators for each mode 
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as 

Xj = a d + a], Y d = -i(aj - a]), (1) 
so that the Heisenberg uncertainty principle requires V(Xj)V(Yj) > 1. 

A. Three-mode Einstein-Podolsky-Rosen correlations 

The EPR argument was introduced in 1935 in an attempt to show that quantum mechan- 



ics could not be both complete and consistent with local realism 







Schrodinger replied 



that same year by introducing the concept of entangled states which were not compatible 



with c 
mond 



assical notions such as local realism 



3id H 



121. In 1989 Reid 1 131, and Reid and Drum- 



1^ | proposed a physical test of the EPR paradox using optical quadrature amplitudes, 
which are mathematically identical to the position and momentum originally considered by 
EPR. Reid later expanded on this work, demonstrating that the satisfaction of the 1989 
two- mode EPR criterion always implies bipartite quantum entanglement Q]. Tan made a 
similar demonstration in the context of teleportation, considering the outputs from a non- 
degenerate optical parameteric amplifier (OPA) mixed on a beamsplitter |l6j . In this article 
we use an extension of Reid's original approach to the case of tripartite correlations, where 
quadratures of three different optical modes are involved. This extension was developed and 
formally proven to demonstrate the presence of tripartite entanglement by Olsen, Bradley 
and Reid in Ref. Q], so that we shall call these the OBR criteria. 

There are two ways to consider the experimentally accessible form of the OBR criteria, 
depending on whether we use information from two quadratures to infer properties of the 
other, or information from one to infer combined properties of the other two. In the first case 
we make a linear estimate of the quadrature Xi from the properties of the combined mode j+ 
k, using parameters which can be optimised, both experimentally and theoretically jl3. \v\. 



It has been shown jl OL Il4j that minimising the root mean square error in this estimate leads 
to an optimal inferred variance, 

v(x,x j ±x k )} 2 

V mf (Xi) = V(Xi) - ± s s i -, (2) 

V(X j ± X k ) 

where V(A, B) = (AB) — (A)(B). We follow the same procedure for the Y quadratures to 
give expressions which may be obtained by swapping each X for a Y in the above to give 
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the optimal inferred estimate 



V(Yi,Yj±Y k 

^m-vm- v(Yj±Yk) 

A demonstration of the EPR paradox can be claimed whenever theory predicts 



(3) 



V inf (Xi)V inf (Yi) < 1. (4) 

As was proven jioj], this demonstration for the 3 possible values of i is then sufficient to 
establish tripartite entanglement, without any assumptions having been made about whether 
the states involved are Gaussian or not. 

Following the same logic, if we use the properties of mode i to infer properties of the 
combined mode j + k, we find that there is a demonstration of the other three mode form 
of the EPR paradox whenever 



V mf (X 3 ± X k )V inf (Yj ± Y k ) < 4, 



(5) 



where, for example, 



V inf (X,±X k ) = V(X J ±X k 



V{Xi,Xj±X k 



V(X t 



(6) 



As above, this demonstration for the 3 possible combinations also serves to establish complete 
inseparability of the density matrix. 



B. The van Loock-Furusawa inequalities 

A set of conditions which are sufficient to demonstrate tripartite entanglement for any 
quantum state have been derived by van Loock and Furusawa [9|. Using our quadrature 
definitions, the van Loock-Furusawa conditions give a set of inequalities, which we shall refer 
to as the VLF inequalities, 



Via 


= V(X 1 


-X 2 ) + V(Y 1 + Y 2 + g 3 Y 3 ) 


>4, 


Via 


= vix, 


- X 3 ) + V(Y + g 2 % + %) 


>4, 


V 23 


= V(X 2 


-X^ + V^Y + Y + Y) 


>4, 
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where V(A) = (A 2 ) - (A) 2 and the 9i are arbitrary real numbers. As shown in reference Q, 
the violation of the first inequality still leaves the possibility that mode 3 could be separated 
from modes 1 and 2, but this possibility is negated by violation of the second. Therefore, 
if any two of these inequalities are violated, the system is fully inseparable and genuine 
tripartite entanglement is guaranteed. We note also that genuine tripartite entanglement 
may still be possible when none of these inequalities is violated. 

We will now investigate optimisation of the VLF criteria, using the freedom allowed in 
the choice of the gi, which are arbitrary real parameters. A simple minimisation of the 
right-hand sides of Eq. with respect to the gi gives 

-{(%%) + {%%)) 



9i 



92 



■((¥,%) + {%%)) 
(Y 2 2 ) 



93 = — W) — • (8) 

The required variances can now be written as, for example, 
V(X\ - X 2 ) = (X 2 ) + (X 2 ) - 2(X 1 X 2 ), 
V{Y X + % + g 3 Y 3 ) = (Y 2 ) + (Y 2 2 ) + ^(F 3 2 ) + 2 [(Y&) + g 3 ((f^) + (%%))] ■ (9) 
Once this optimisation process has taken place, we find that, for example, 

ViX^Y + Y)} 2 

V(Y 1 + Y 2 + g 3 Y 3 ) = V(Y 1 + Y 2 ) - ± j -, (10) 



where we recognise the right hand side as an inferred variance as introduced in Ref. [10J 
to demonstrate the EPR paradox for three modes, and referred to above. The optimised 
correlations can now be written as 

V 12 = V{X 1 - X 2 ) + V inf (Y + Y 2 ) > 4, 
V 13 = V(X X - X 3 ) + V mf (Y + Y 3 ) > 4, 

V 23 = V(X 2 - X 3 ) + V m f{% + %) > 4. (11) 

We see that the VLF criteria now have the same form as the Duan and Simon criteria for 
bipartite entanglement Q. but with the actual variance V(Yj + Y) replaced by the 
inferred variance V m f(Yj + Y^) of Eq. E| We note that the violation of two out of three of 
the inequalities is sufficient to demonstrate full inseparability. 



III. SYSTEM AND EQUATIONS OF MOTION 

The interaction Hamiltonian used by Bondani et al. jjij uses an undepleted pumps ap- 
proximation and we will begin with a more complete form which quantises all the interacting 
fields. The Hamiltonian describes the coupling of five modes of the electromagnetic field 
in a phase-matched simultaneous sum frequency generation and downversion process in a 
manner analogous to the schemes considered by Olsen and Bradley |7[, and has previously 
been investigated by Ferraro et al. [5( and Smithers and Lu 

Q- In Ref. Q, this five-mode 

Hamiltonian is written as 

Hint = ih (xAa{al + X2a 5 aia 3 } + h.c, (12) 

once we change the indices to agree with Bondani et al. jl| and set the coupling coefficients 
to be real. Due to energy conservation, lo\ = ui\ + U3 and o> 2 = cj 3 + u^, and the necessary 
phae-matching conditions are covered in Ref. Q]. This Hamiltonian approximately describes 
a downconversion process cascaded with a sum-frequency generation process where one of 
the downconverted modes becomes an auxiliary pump mode for the frequency generation 
process. It gives a simplified description because it does not include effects such as dispersion 
within the nonlinear medium, for example. A more accurate method of analysing these types 
of processes has been given by Raymer et al. [21], but the approximations we are using do 
serve to set upper limits on the s quee zing and entanglement available from a more realistic 



treatment of the physical process 



221. 



However, given the above caveat, it is instructive to examine the analytical solutions 
which may be obtained using an undepleted pumps approximation as without a cavity the 
interaction strengths tend to be small and this approximation is generally very accurate. 
Setting K\ = Xi(a4(0)) and k 2 = X2(o5(0)) as real positive constants, the Hamiltonian may 
be written 



'Hint. — ih 



K\ fojaj — dih^ + K 2 ^2^3 — ^2^3^ j (13) 



from which we find the Heisenberg equations of motion, 

ddi i 

IT = 

d(l2 

-77- = K 2%, 

dt 



For later convenience we will rewrite the above as equations of motion for the quadrature 
operators, finding 

dX x 



dt 
dYx 
~dt 
dX 2 

dt 
dn 

dt 
dX 3 

dt 
dYz 

dt 



-K1Y3, 

K 2 X 3 , 

K±Xi — K2X2, 

- Kl Yi - k 2 Y 2 . (15) 



These equations can now be solved analytically to give the solutions for the operators as 
functions of their initial values, which will all be zero for this system. However, due to 
bosonic commutation relations, not all the moments vanish and at t = with all the output 
fields as vacuum, we have (Xi(0)Xj(Q)) = (Yi(0)Yj(0)) = 5y. This is all the information we 
need to find useful time-dependent solutions for the variances and covariances needed for 
the correlations which establish tripartite entanglement. 



A. Hyperbolic solutions 

We find that there are three classes of solutions for different regimes, depending on 
whether k\ > n\, k\ < k\ or k\ = k\. The last of these was treated in Ref. \l\ and we will 
not consider it further here. For k\ > k x , Q is imaginary and the solutions are periodic, 
while for k 2 < k\ the solutions are hyperbolic. We will begin with the correlations for the 
hyperbolic solutions, as this is the operating regime of the Bondani experiment 

Setting = a//c? — k 2 , we find these solutions as 

Mt) = "1 cosh ^ ~ ^ (0) - ^ (co g"* ~ ^ (o) + ^g^x 3 (o), 

~ k\ cosh Qt - kI * K1K2 (cosh Qt - 1) y Atx sinh Q.t y 

' ' ~ 9? ^ ' Q? Y2 ^> Q y 3lUJ, 



-v . . Ki sinh fit - . , k 2 sinh fit - , s ^ . . , ^ 
X 3 {t) = Xi(0)--?-^ X 2 (0)+X 3 (0) cosh fit, 

- . . Ki sinh fit - . . Ko sinh fit - . , - . . , „ . 
Y 3 (t) = -^—^ Fx(0) - y 2 (0) + y 3 (0) cosh fit, (16) 

which contain all the information needed to calculate the VLF and OBR correlations in 
the approximations we are using, except in the case where In this case the 

above solutions are not well defined but the equations may still be solved using stochastic 
integration, as was done in Ref. For k\ > k\, the time-dependent moments which we 
need are 

2k 2 

(Xl) = (Y?) = 1 + [k\ sinh 2 fit + 2k 2 (1 - cosh fit 
(If) = (¥?) = 1 + ij [2^ 2 (coshfit - l) 2 ] 

(XiX 2 ) = -(FiV2> = [(«? + «2) ( cosh 1) 2 + ^ 2 sinh 2 fit] 



fi 4 

(XxXs) = -(Y 1 Y 3 } = 2KlS ^ aQ \ K lcoshnt-4) 



(X 2 X 3 ) = (Y 2 %) = (cosh fit - 1) sinh fit. (17) 

We note here that the above expectation values are actually the variances and covariances 
for single modes, as the expectation values of the amplitudes are all zero. 



B. Periodic solutions 



We will now look at the case where k 2 > k 2 , which leads to solutions expressed in terms 



of periodic functions. Setting £ = \J k 2 — k 2 , we find 

Mt) = Kl - f^^ xm + ^ 2(c y ~ ^ (o) + Kl 

s s s 



A = KlK2 (i - c°sa) li(0) + - 4 Mo) + *gg&( ), 



/ti/t 2 (cos£t- 1) ~ K 2 cos£t- k 2 - K 2 sin£t, 

£2" y H U ^ + £2 y 2V U J + £ 



- , , Ki sin £t - , , k 2 sin £t . . - . . 
X 3 (t) = - q Xi(0) - - q X 2 (0) + X 3 (0) cos£t, 

- . , k,\ sin £t - . . /c 2 sin ft - . . - . . . 
Y 3 (t) = - Ki(Q) - ^ y 2 (0) + y 3 (0) cos£t, (18) 
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FIG. 1: (Colour online) The analytical solutions of the VLF correlations, with k\ = 1.2/^2- Any two 
of the correlations falling below 4 is sufficient to demonstrate that genuine tripartite entanglement is 
present. The solid lines use the optimised expressions. All quantities shown in these and subsequent 
graphs are dimensionless. 

which lead to the solutions for the moments, 

,* 9% 2k\ \2k\(\ -cos it) - k\ sin 2 ft] 
(X 2 ) = (>f> = 1 + U ^ l - il 

, a,, 2ft 2 K 2 ,(cOs££ — l) 2 

{X X X 2 ) = -(%%) = ^ [{k\ + - cos ^ - K 2 sin 2 £t] , 

(XiXa) = -(y^) = p [2« 2 sm£t-/c 2 sin2£t] , 



, - - , 2/t? Ko sin ft r 
(X 2 X 3 ) = (Y 2 Y 3 ) = 1 2 3 ^ [1 - cos£t] • (19) 



IV. ENTANGLEMENT RESULTS 

Analytical expressions can be found for both the VLF and OBR correlations using the 
results of Eqs. [T7| and HH1 but as these can be rather unwieldy we will present our results 
graphically. In the interests of compact notation we will use the shorthand 14,- for the corre- 
lation which contains V(Xi — Xj). In Fig. ^we show the results of the VLF correlations for 
the hyperbolic solutions, with K\ = 1.2re 2 - The dash-dotted lines are the basic expressions, 
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FIG. 2: (Colour online) The analytical solutions of the VLF correlations, with k 2 = 1.8«i. Any two 
of the correlations falling below 4 is sufficient to demonstrate that genuine tripartite entanglement 
is present. The solid lines use the optimised expressions. 

without any optimisation, and demonstrate that genuine tripartite entanglement is present 
over a small range of interaction strength. The solid lines (of the same colour online) are the 
expressions optimised as in Eq. ^Jand are seen to violate the inequalities over a wider range. 
Perhaps the effect of this optimisation is that it allows for the demonstration of entanglement 
as soon as the interaction is non-zero, whereas the expressions without optimisation need 
some finite interaction before any of them go below 4. This is not a contradiction as entan- 
glement may be present even if the inequalities are not violated, in contrast to the Duan 
and Simon criteria for Gaussian bipartite systems, which provide necessary and sufficient 
conditions jlfl . We are not aware of any criteria for tripartite continuous- variable entan- 
glement which provide both necessary and sufficient conditions. In Fig. El we present results 
for the same correlations in the regime of periodic solutions, with k 2 = 1.8«x- We again 
see that, as expected, the optimisation procedure allows for demonstration of entanglement 
over a wider range of interaction strengths. 
The three-mode EPR correlations, 

OBR, = V inf (Xi)V inf (Yi), 
OBRj k = V^'iXj + X^V^iYj + Y k ), (20) 
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FIG. 3: (Colour online) The analytical solutions of the OBR two-mode inference correlations. 
On the left hand side, K\ = 1.2^2, while on the right hand side K2 = 1.8«i. Although all three 
correlations should be below one to demonstrate genuine tripartite entanglement, in each case only 
one of them (OBRi) goes below this level, while the other two are both exactly equal to one. 

may be expressed in terms of the operator moment expectation values using 

[(XiXj) + (XX k )f 



V mf {X) = {X?) 



V inf {%) 
V mf {X 3 +X k ) 



(m - 



(Xf) + (Xi)+2(X j X k y 
[(YjYj) + (Y&)] 2 



(X") + (Xl) + 2(X j X k ) - 



[(XX,) + {XiX k )Y 



V^(Y 3 + Y k ) = (Yf) + (Y?) + 2(Y j Y k ) - 



(X?) 
[(YjYj) + 

(V) 



(21) 



In Fig. El we give the results for the OBRi in both the periodic and hyperbolic regimes. 
Neither of these results, which come from inferring the properties of a single quadrature 
from the properties of a combined two-mode quadrature, gives evidence of genuine tripartite 
entanglement. In fact, all that these particular correlations succeed in demonstrating is 
that the combined density matrix, p 123 , cannot be separated in the manner p 12 3 = P1P23, 
while leaving open the possibilities P123 = P2P13 and P123 = p-sPu- This shows that choosing 
to measure these particular criteria to demonstrate entanglement would not be sensible, in 
contrast to the triply nonlinear system considered in Bradley et al. [f| and Olsen et al. 



e, m 
[10]. 
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where the symmetries of the interaction Hamiltonian meant that any choice of the VLF or 
OBR criteria was equally useful, with all three giving comparable results. 

However, we do find that with the present system the three-mode EPR correlations 
(OBRij), which are defined using the properties of one quadrature to infer properties of 
a combined quadrature which involves the other two modes, are operationally useful. As 
shown in Fig. H] and Fig. there is an unambiguous demonstration of the inseparability of 
the density matrix almost as soon as the interaction begins. This demonstration continues 
well past the point where the undepleted pumps approximation is expected to lose its va- 
lidity. Hence, if entanglement were to be demonstrated experimentally with this scheme, 
measurement of these three correlations would be the preferred option. 




0.2 0.4 0.6 0.8 1 1.2 



FIG. 4: (Colour online) The analytical solutions of the OBR correlations which infer combined 
mode properties from those of a single mode, with k\ = 1.2^2- All three correlations should be 
below four to demonstrate genuine tripartite entanglement. 
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FIG. 5: (Colour online) The analytical solutions of the OBR correlations which infer combined 
mode properties from those of a single mode, with K2 = 1.8k±. All three correlations should be 
below four to demonstrate genuine tripartite entanglement. 

V. CONCLUSIONS 

We have examined the interlinked interaction scheme of Bondani et al. |l] in terms 
of its suitability for producing output fields which exhibit genuine tripartite entanglement. 
Using the undepleted pumps approximation, which is valid for small interaction strengths, we 
have calculated correlations using three different approaches, one which uses the properties 
of combined quadratures and may be optimised, and the other two which use three-mode 
generalisations of the EPR argument. We find that these correlations give different answers 
to the question of whether tripartite entanglement is present in a particular regime and 
that the most sensitive is that developed by Olsen et al. Q| to infer the properties of a 
combined mode from those of a single mode. These inequalities are successful in detecting 
the entanglement over a large regime where the VLF criteria give a false negative. The 
fact that the correlations give different answers is not contradictory as they all provide 
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sufficient but not necessary criteria and this is a good example of how investigations of 
continuous variable entanglement become more complicated once we have more than two 
modes involved. 
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